Gel to glass transition in simulation of a valence-limited colloidal system 
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We numerically study a simple model for thermo-reversible colloidal gelation in which particles 
can form reversible bonds with a predefined maximum number of neighbors. We focus on three 
and four maximally coordinated particles, since in these two cases the low valency makes it possible 
to probe, in equilibrium, slow dynamics down to very low temperatures T. By studying a large 
region of T and packing fraction </!> we are able to estimate both the location of the liquid-gas phase 
separation spinodal and the locus of dynamic arrest, where the system is trapped in a disordered 
non-ergodic state. We find that there are two distinct arrest lines for the system: a glass line at high 
packing fraction, and a gel line at low (j> and T. The former is rather vertical ((^-controlled), while the 
latter is rather horizontal (T-controUed) in the {(p — T) plane. Dynamics on approaching the glass 
line along isotherms exhibit a power-law dependence on 0, while dynamics along isochores follow 
an activated (Arrhenius) dependence. The gel has clearly distinct properties from those of both a 
repulsive and an attractive glass. A gel to glass crossover occurs in a fairly narrow range in (p along 
low T isotherms, seen most strikingly in the behavior of the non-ergodicity factor. Interestingly, 
we detect the presence of anomalous dynamics, such as subdiffusive behavior for the mean squared 
displacement and logarithmic decay for the density correlation functions in the region where the gel 
dynamics interferes with the glass dynamics. 

PACS numbers; 82.70.Gg, 82.70.Dd, 64.70.Pf 



I. INTRODUCTION 

Systems composed of mesoscopic solid particles dis- 
persed in a fluid are named colloids. The properties 
of the particles and of the fluid can be controlled via 
chemical or physical manipulations to a great extent. 
As a result the particle-particle interaction can be tuned 
from very short range depletion attractions to very long 
range Coulombic repulsion P|, making colloids impor- 
tant both in terms of basic scientific research and in- 
dustrial applications |E IE ^ IS- The canonical model 
system for colloids is the hard sphere system, for which 
sterically stabilized colloidal particles such as PMMA 
provide a very accurate experimental realization |^, Q ■ 
Hard sphere colloids have been used to directly observe 
crystal nucleation and to test theoretical predictions sur- 
rounding glass transitions The addition of small, non- 
adsorbing polymers to a hard-sphere solution leads to a 
short-range effective attraction between colloids through 
the so-called depletion interaction 0]. The size of 
the small polymer controls the range of attraction, while 
the concentration controls the strength of attraction uq. 
Neglecting the role of solvent interactions, these sys- 
tems can be simulated on a computer with a short-range 
attractive potential, as simple as a hard core comple- 
mented by a square well (SW). Colloid-polymer mix- 
tures have been found to offer new scenarios of arrested 
states. These hard sphere plus short range attraction 



systems exhibit the usual hard sphere glassy dynam- 
ics near (j) ~ 0-6. However, a fascinating phenomenon 
arises when the range of attraction is less than approx- 
imately 10% of the hard sphere diameter, that is a re- 
entrance of the glass transition line, predicted by Mode 
Coupling Theory {MCT)_ElElEkiand confirmed by 
several experime nts MSl lid llTl llSl and simula- 
tions [13, In m lillllPin this case, for a particular 
range of (f>, arrest can be achieved by either increasing or 
decreasing T/uq, the ratio of temperature T to attraction 
strength uq. The two types of glass are now commonly 
named hard-sphere or repulsive glass for the one at higher 
T, and attractive glass for the one at lower T. 

At low densities short-range attractive colloids exhibit 
particle clustering and gelation [25l Wll [2^ . Recently, 
several numerical works have focused on colloidal gelation 
[23|33,mi32„,33^^J, with the aim of better characteriz- 
ing colloidal gels and attemptingto formally connect gel 
to glass formation [H El IM iM HI ■ Following MCT 
ideas, the gel state was interpreted as a low (p extension of 
the high cj) attractive glass ^jHHHIs^l and, generating 
some confusion, the term " gel" is still often interchanged 
with the term 'attractive glass'. 

Differently from chemical gelation , th at was exten- 
sively studied in polymer phvsics's S 1391 140| and mod- 
eled in computer simulations 41. 42. i43l |44L l45l 11^ . col- 
loidal gelation is still quite poorly understood. Gelation 
arises when a stable particle network is formed due to 
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bonding. For chemical gels, bonds are irreversible, and 
thus gelation can be explained in terms of percolation 
theory. However, the bonds intervening in colloidal ag- 
gregation have an energy typically of the order of fesT, 
as for example bonds induced by depletion interactions. 
Thus such bonds are mostly transient at finite attrac- 
tion strength. The existence of a finite bond lifetime 
creates a gap between the location of the percolation line 
and the dynamic arrest (gel) line. Earlier simulations 
on a lattice]^ 0, 01 and ofi^-lattice^^ have discussed 
the influence of reversible bond formation on the gelation 
process. 

The quest for bond stabilization often calls for explor- 
ing regions of very high attraction strength in the phase 
diagram to increase the bond lifetime and promote gela- 
tion. However, in this region, at low packing fraction 
(j) and at low T (or at high uq values), a phase sep- 
aration into gas (colloid-poor) and liquid (colloid-rich) 
always takes place|3|- Recent theoretical studies have 
addressed the question of the relative location of the 
phase separation and of the attractive glass line. It has 
been found J£, j^Qj that the attractive glass line meets 
the spinodal at a finite temperature, on the high density 
side. At low </>, arrest in (spherically interacting) short- 
range attractive colloids occurs only as arrested phase 
separation, where the liquid phase is glassv, and hence 
the phase separation cannot proceed fully |43 ■ This sce- 
nario was found in numerical simulations of the square 
well potential for a well width ranging from 15% to ar- 
bitrarily small values Hfj i5^1l| . as well as in Lennard- 
Jones potentialjs^- Gels as a result of interrupted phase 
separation have been identified in experiments 53] and 
simulations [49II53 . l55i | . For extremely deep quenches (at 
very low T) irreversible gelation rnay occur through dif- 
fusion limited cluster aggregation [Sg, [S^, IH, |5^ |60| . 

Other mechanisms must be invoked for stabilization 
of thermo-reversible gels and suppression of the spinodal 
line'oi'l. A possibility is to consider the effect of resid- 
ual charges on the colloidal particles, that may produce 
a longe-range repulsive barrier in the interparticle poten- 
tial, thus efficiently stabilizing the bonds and preventing 
the condensation of the liquid phase. The emergence of 
an equilibrium cluster phase has been recently evidenced 
both in experiments |6al63| and simulations |30l l3lL l3a | on 
charged colloidal and protein suspensions. A dynamical 
arrest transition then follows, driven by electrostatic re- 
pulsions between the clusters ,30,] or by cluster branching 
and percolation |^ 1^3 • 

Very recently, forming a gel in equilibrium, without 
the interference of phase separation, has been achieved by 
limiting the number of attractive interactions (bonds) be- 
tween nearest neighbors. Following ideas first introduced 
by Speedy and Debenedetti 64, 65], a model for thermo- 
reversible gelation has been studied numerically^ . This 
model, in which particles interact via a SW potential 
with the addition of a geometrical constraint in the max- 
imum number of bonds rimax a particle can form, de- 
scribes particles interacting via a hard-core plus n^ax 



randomly- located sticky spots [63 ■ In this respect, the 
model retains the spherical symmetry but incorporates 
features of associating liquids [6^. Other possibilities 
to limit the number of bonded pairs implementing ge- 
ometrical constrains and retain ing a spherical pair po- 
tential can be found in Ref. [681 l69j , while the same ba- 
sic ideas have inspired recent work on self-assembly of 
super- molecular/nano structures [toI Flf . In the liter- 
ature it is possible to find several related models, also 
based on limited short-range attractions in which the 
bonding constraint is imposed via angular degrees of free- 
dom, though the focus has not always been on their dy- 
namic properties [H ItIItIItI Fsl] . In general, ref | 33| 
puts forward the hypothesis that, only when a restricted 
part of the colloidal surface is active in the formation of 
attractive bonds, dynamic arrest at low (f> can be observed 
in the absence of a phase separation. An experimental 
test of the previous hypothesis will hopefully be provided 
by the new generation of 'patchy' colloids, or colloids with 
'sticky spots' [76|, which is about to be synthesized. 

Understanding gelation at low densities in short-range 
attractive models may also be relevant to the study of 
proteins, since they are expected to belong to the class 
of short-range attractive interacting systems|23,l23- In- 
deed, the formation of arrested disordered states at low 
densities often interferes with crystallization, and this is 
possibly one of the reasons why proteins are often diffi- 
cuh to crystalhze [T^Fflj ]. 

In a recent Letter '33l] , we showed that the nmax model 
(detailed in section II below) allows us to study thermo- 
reversible gels. We showed that the signatures of the gel 
state, as for example the non-crgodicity factor /q, are 
quite distinct from those of both the hard sphere and the 
attractive glass. In the present study, we explore the full- 
ed dependence of dynamics at low T for both rimax = 3 
and Umax ~ 4. We observe a transition over a small 
range in (f) from a gel to a repulsive glass. In Section 
II, we give details of the model and simulations. Section 
III contains the results for the calculated phase diagram 
and compares the relative location of the thermodynamic 
and of the kinetic arrest lines. We also report static and 
dynamic correlation functions. In Section IV we discuss 
results and in Section V conclusions are drawn. 



II. SIMULATION DETAILS 

We perform event-driven molecular dynamics simula- 
tions of iV = 10000 particles of mass m = 1 with di- 
ameter a — 1 (setting the unit of length) interacting 
via a limited-valency square well potential. The depth 
of the well uq is fixed to 1, and the width A of the 
square well attraction is such that A/(cr -f A) = 0.03. 
T is measured in units of uq, and the unit of time t is 
(7{m/uoy^'^ ■ This system is a one-component version of 
the binary mixture that has been extensively studied pre- 
viously [23, Isol Isil ] . In the following we will use the 
acronym SW to indicate the A/ (cr -I- A) = 0.03 standard 
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square well potential. The limited-valency condition is 
imposed by adding a bonding constraint. The square 
well form can be used to unambiguously define bonded 
particles, i.e., particles with centers lying within a and 
cr + A of each other are bonded. The interaction between 
two particles each having less than nmax bonds to 
other particles, or between two particles already bonded 
to each other, is thus given by a square- well potential. 



oo r < cr 
r>a + A. 



(1) 



When i and/or j are already bonded to ri„iax neighbors, 
then Vij{r) is simply a hard sphere (HS) interaction. 



oo 




r < a 
r > a. 



(2) 



The resulting Hamiltonian of the system has a many- 
body term containing information of the existing bonds. 
Due to the fact that the list of existing bonds is necessary 
at any instant of the simulation, configurations are saved 
storing also the bond list. Moreover, all simulations are 
started from high temperature configurations where all 
particle overlaps are excluded within the attractive well 
distance cr -I- A. In cases of multiplicity of possible bond- 
ings, such as for example when a bond is broken for a 
particle that was fulfilling the nmax allowed bonds and 
more than one neighbour particle lie within its attractive 
well, a random neighbour, with less than n^ax bonds, is 
chosen to form the new bond. 

The idea of constraining the number of square well 
bonds a particle can form was introduced by Speedy and 
Debenedetti 0, |6^ . In contrast to their original version 
of the model, where triangular closed loops were not al- 
lowed, in the present model no constraints on minimal 
bonded loops are introduced. Our model can be con- 
sidered as a realization of particles with nmax randomly 
located 'sticky spots' [g^. The Hmited- valency properly 
defines the ground state of the systems, corresponding 
to a potential energy per particle — uo'^max/2. This is 
achieved when every particle has rimax filled bonds. We 
note that fully bonded clusters of finite size may occur. 
The smallest fully connected cluster sizes are: 4 (tetra- 
hedra) for rimax = 3; 6 (octahedra) for rimax = 4; and 12 
(icosohedra) for rtmax = 5. This introduces the intriguing 
possibility of forming a hard sphere gas of such clusters 
at low (j) and T, as seen from the simulations. We study 
in depth the cases nmax = 3 and 4, for which we already 
know that there exist significant differences in the loca- 
tion of the spinodal lines as compared to the rimax = 12 
(or standard SW) case. 

For all state points simulated, we first equilibrate the 
system at constant T until the potential energy and pres- 
sure of the system reach a steady state, and the MSD 
reaches diffusive behavior, i.e., during the equilibration 



time particles move on average at least one particle di- 
ameter. A subsequent constant energy simulation is used 
to gather statistics for all reported quantities. 

To estimate the equilibrium phase diagram, we cal- 
culate the gas-liquid spinodal and the static percolation 
line. The latter is defined as the locus in ((/>, T) such that 
50% of the configurations possess a spanning, or perco- 
lating, cluster of bonded particles. To characterize the 
structure and dynamics of the system, we evaluate: the 
static structure factor, 



(3) 



Siq)^(\pM\ 



the mean squared displacement (MSD), 



N 



< r^(t) >^< ^ \m - /N >, (4) 



the diffusion coefficient. 



D = Yim ( - I^Nt, (5) 



i 1=1 



the dynamic structure factor, or density autocorrelation 
function. 



i^,(t)^(p,(iV_,(o))/(p,(o)p_,(o)) 



(6) 



and its long time limit or plateau value fq, i.e. the non- 
ergodicity parameter. In all cases, (.) denotes an ensem- 
ble average, ri is the position vector of a particle, q is 
a wavevector and i labels the N particles of the system, 
while pq{t) = J2fLi exp {-iq- n)- 

Also, we monitor the bond lifetime correlation func- 
tion, averaged over different starting times, and defined 
as, 

Mt) ^ {Y.n,jit)n,,{0))/[NBiO)l (7) 

i<j 

where (t) is 1 if two particles are bonded up to time 
t and otherwise, while Nb{0) = {^i^j nij{0)) is the 
number of bonds at t = 0. We note that (f>B counts 
which fraction of bonds found at time t = persists 
at time t, without ever breaking within the store rate 
of configurations. Associated with c/fsit), we extract an 
estimate of the bond lifetime tb via stretched exponential 
fits. 

Although colloidal systems are more properly modeled 
using Brownian dynamics, we use event driven molecu- 
lar dynamics due to its efhciency in the case of step- wise 
potentials. While the short-time dynamics is strongly af- 
fected by the choice of the microscopic dynamics, the long 
term structural phenomena, in particular close to dynam- 
ical arrest, are rather insensitive to the microscopic dy- 
namics [73, |82| . To have a confirmation of this, we also 
performed additional simulations where the effect of the 
solvent was mimicked by so-called ghost particles 's^. 
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In particular, we studied a system of 1000 colloidal par- 
ticles and 10000 ghost particles and we found that the 
long-time behavior of the dynamical quantities, such as 
the T dependence of tb and the q dependence of fq are 
independent from the microscopic dynamics. We also 
note that equivalence between Newtonian and Brownian 
dynamics is not guaranteed when studying small length 
scales (as for example the decay of density fluctuations at 
large 9), i.e. for intra cage motion, since there the micro- 
scopic dynamics may affect the shape of the decay of the 
correlation function. This is most relevant in gel systems 
for which the cage length can be significantly larger than 
the particle size. 



III. RESULTS 

In this section we examine the results of our simulation 
in terms of the thermodynamic and dynamic quantities 
mentioned above. We focus our attention on the cases 
'T-max = 3 and rimax = 4, where a significant suppression 
of the li quid -gas spinodal as compared to the SW case is 
observed|33j. The suppression of the critical temperature 
and the shrinking of the unstable region in the {(j) — T) 
plane makes it possible to study in one-phase conditions 
state points characterized by an extremely slow dynam- 
ics. More precisely, it is possible to study without en- 
countering phase separation all (p > 0.2 for rimax = 3 and 
(j) ^ 0.30 for rimax = 4. In the SW case, phase separa- 
tion was encountered already at T « 0.32, when particle 
mobility is always large 0. 

A. Phase Diagram 

Ref. la E3 

show that constraining the number of 
bonded neighbors reduces the energetic driving force for 
particle clustering. Therefore, as nmax decreases, the 
phase separation transition (that we monitor by studying 
its spinodal) , and the percolation line both shift to lower 
T (along isochores) and lower (f) (along isotherms). 

Fig. ^ reports the percolation and spinodal loci for 
JT-max = 3 (a) and n,„ax = 4 (b). To evaluate the lo- 
cation of the spinodal line we interpolate the pressure 
P{4)) and search for the condition dP/dcj) — 0. To bet- 
ter track down the spinodal and estimate the effect of 
the spinodal in its vicinity we also determine the loci of 
constant S{q — > 0), that we name 'iso-S'(O) lines', that 
can be considered as precursors of the spinodal line. In- 
deed, S{Q) is connected to the isothermal compressibil- 
ity KT = {d(j)/dP)/(t) by 5(0) = pkeTnT. The iso-S'(O) 
lines are calculated as loci of constant {dP/d(p)/kBT. We 
cross-check these estimates with the less precise value ob- 
tained calculating directly S{q 0). The iso-iS'(O) lines 
are shown to emphasize that, in bonded systems, an in- 
crease at small-q in the scattering intensity can arise in 
the one-phase region due to the vicinity of the spinodal 
curve. 




Ql , 1 , 1 , 1 , 1 , 1 , 1 

0.1 0.2 0.3 0.4 0.5 

FIG. 1: Phase diagram for rimax = 3 (a) and rimax = 4 
(b), showing the spinodal (dashed lines with squares), per- 
colation (solid lines with open circles), iso-diffusivity loci 
where DjyfT — constant (lines with triangles), 'iso-S'(O) 
lines' (dashed lines). Also shown is the extrapolated glass 
line, labeled as (pc from power-law fits, see text. 



Fig. n] also shows iso-diffusivity lines, i.e. lines where 
D/^/T is constant. The scaling factor y/T is used to take 
into account the trivial contribution of the thermal veloc- 
ity with T. The investigated values of D/^/T cover four 
orders of magnitude. Such iso-diffusivity lines are pre- 
cursors of the dynamical arrest transition, corresponding 
to Z3 ^ 0. Previous works pH l2^ Is^] have shown that 
these lines provide estimates of the shape and location 
of the arrest line. The isodiffusivity lines show an in- 
teresting behavior. They start from the high-0 side of 
the spinodal curve and then end up tracking the high-T 
hard-sphere limit. They are rather horizontal (parallel 
to the (/)-axis) at low (j) rapidly cross to a vertical 
shape (parallel to the T-axis) at high (j). The crossing 
from horizontal to vertical becomes sharper and sharper 
on decreasing D/^/T. In the SW case, the iso-diffusivity 
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FIG. 2: Power law fit of inverse difTusivity D = ^((^ic — 
<j>) for (a) rimax = 3 and (b) nmax = 4. Lines are guides 

to tlie eye. Insets (a-i) and (a-ii) sliow tiie beliavior of (pc in 
the (j> — T plane. Inset (b-ii) and (b-ii) show the exponent 7 
as a function of T. 



lines exhibit a reentrance in 0, in agreement with the 
predictions of the Mode Couphng Theory (MCT). The 
reentrance becomes more and more pronounced at lower 
and lower D/^/T values. In the rimax case, a reentrant 
shape is hardly observed. Indeed, in the SW the reen- 
trance arises from the competition between cages created 
by the nearest neighbors excluded volume, with a typical 
hard-sphere localization length ^ O.Ict, and cages created 
by bonding with a localization length ^ A. In the rtmax 
case, such a competition becomes less effective due to the 
smaller number of bonds. 



B. Estimation of dynamical arrest lines 

To provide an estimate of the dynamical arrest lines, 
we can identify a range of parameters where the charac- 
teristic time follow a power-law dependence. In the case 
D is selected, data can be fitted according to D~^{<j)) — 
A{cj)c — <jji)~''^'^\ where 0c is the best estimate for the glass 
line. Fig|21shows in log- log scale D^^{(f>) / ^/T vs {(f>c — 0) 
for riniax = 3 and rimax — 4, where the 0c values are cho- 
sen by a best-fit procedure. For (0c — 0) ^ 0.3, data are 
found to be well-represented by power-laws for at most 
two decades in D. The fits are performed over the range 



(0c — 0) < 0.3. The T dependence of the fit parameters 
is also reported in Fig[21 The two fit parameters vary 
almost in phase with each other for both nmax values. 
At high temperatures {T > 0.3) attraction does not play 
a role, the arrest line is almost vertical and 0c and 7 are 
practically constant. For T < 0.3 smaller values of 0c and 
7 are found. The values of 7 are rather small at high T, 
close to the lowest possible value allowed by MCT, and 
become smaller than the lowest possible value allowed 
by MCT at low T. We note on passing that crystal- 
lization limits the range over which dynamic measure- 
ments in (metastable) 'equilibrium' can be performed. 
The region where we detect crystallization varies with T. 
At high T crystallization happens already for > 0.54, 
while at low T, crystallization does not intervene up to 
= 0.56. Thus, we cannot fully rely on the high T fits 
as (0 — 0c) is always large. Notwithstanding this, the 
estimate 0c ~ 0.58 is reasonable. On the other hand, at 
low T, the vicinity to 0c increases, but the range of D 
values over which a power-law can be fitted decreases to 
about one decade, making the fit questionable. The fact 
that 7 decreases well below the lowest meaningful MCT 
value could tentatively be associated to a difficulty of the 
theory to handle the crossing to an energetic caging (see 
also below). The resulting (j>c{T) line is also reported in 
Fig. n Independently of the fitting procedure, a clear 
vertically shaped arrest line, driven mostly by packing, 
is observed on isothermal compression. 

To provide a better estimate of the arrest line at low 
we have studied the behavior of D with T along iso- 
chores. We have tried two routes. The first one consists 
in performing again power law fits but in temperature 
along an isochore, i.e. D{T) ~ A{T — Tc)"''^ , selecting an 
appropriate T-interval. For both rimax = 3 and 4, such 
fits appear to hold for a rather small interval in D and 
are strongly dependent on the chosen T-interval selected 
in the fit procedure. In this way, we cannot extract a 
reliable estimate for Tc- Indeed, for the same state point 
Tc could vary from 0.2 to 0.08 depending on the fitting 
interval. However, fixing a T-interval of fitting for all 
isochores, the resulting Tc(0) is again rather flat for both 
?^max values. The second route is more robust and is 
based on the observed low T behavior, where data are 
found to follow very closely an Arrhenius law. Fig. O 
shows D as a function of 1/T, for both rimax = 3 and 
4. Arrhenius behavior of D is observed at all at low 
T. The activation energy is around 0.45 for nmax — 3 
and of 0.55 for nmax = 4. Since at the lowest studied 
temperatures the structure of the system is already es- 
sentially T-independent (see later on the discussion con- 
cerning FigO, there is no reason to expect a change in 
the functional law describing the T — > dynamics. In 
this respect, the true arrest of the dynamics is located 
along the T = line, limited at low by the spinodal 
and at high by crossing of the repulsive glass transi- 
tion line. This peculiar behavior is possible only in the 
presence of limited valency, since when such a constraint 
is not present, phase separation preempts the possibility 
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FIG. 3: Arrhenius plot of D for rimax = 3 (a) and 4 (b) along 
all studied isochores without intervening phase separation. 
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FIG. 4: Summary of the thermodynamic and kinetic phase 
diagram for rimax ~ 3, including spinodal (dashed lines with 
filled circles), iso-diffusivity loci where D/VT = 0.005, 0.0005 
(lines with triangles and squares), iso-5'(0) locus where S(0) = 
1.5, extrapolated glass, labeled as <^c, and gel, labeled as D = 
0, lines respectively from power law and Arrhenius fits (Tc = 
0), see text. 



refer the reader to a future work to compare these results 
with corresponding MCT predictions for the same model 



of accessing the T — > Arrhenius window. 

We can summarize the dynamical arrest behavior in 
Fig. ^ One locus of arrest is found at high cj), rather 
vertical and corresponding to the hard-sphere glass tran- 
sition. The isodiffusivity lines suggest a rather flat arrest 
line. Two different loci could be associated to arrest at 
low T. One defined by the Tc of the power-law fits and 
one at T = associated to the vanishing of D according 
to the Arrhenius law. Notwithstanding the problem with 
the power-law fits and the big undeterminacy on Tc, it 
would be tempting to associate the Tc-line to the attrac- 
tive glass line, at least as a continuation of it at low (j), 
and interpret the wide region between the two lines as a 
region of activated bond-breaking processes [s^l- Data re- 
ported in the next sections will show that, in the present 
model, the identification of such line with an attractive 
glass line (or its extension to low density) is not valid, in- 
dependently of the fit results. Indeed, we will show that 
particles are never localized within the attractive well 
width A, at any (f). On the other hand, the establishment 
of a percolating network of long-lived bonds, that we will 
refer to as a gel is identified. Rcf. 33] shows that arrested 
states at low cf) and T are profoundly different from both 
attractive and hard-sphere glasses. Now we are going to 
investigate if on increasing (j) there is a crossover or a 
transition from gel to attractive glass, or whether the at- 
tractive glass exists at all in limited-valency models. We 



C. Static structure factor 

This section reports results for the static structure fac- 
tors for various studied T and and both rimax = 3 and 
'T-max = 4. Results along an isochore and an isotherm are 
general for both studied Jimax values. 

Fig. [SJl shows the evolution with temperature of S{q) 
at the lowest accessible cj) (i.e. the lowest 4> where phase 
separation is not present). On lowering T, S{q) shows 
an increase of the intensity at small wave vectors, which 
saturates to a constant value when most of the bonds 
have been formed. This indicates that the system be- 
comes more and more compressible, with large inhomo- 
geneities, characterizing the equilibrium structure of the 
system. The inhomogeneities can be seen as an echo of 
the nearby phase separation or, equivalently, as a conse- 
quence of building up a fully connected network of par- 
ticles with low coordination number. The large signal at 
small g is a feature of S{q) which is often observed in gel 
samples |23,0|- However, sometimes it may be difficult 
to discriminate between a true equilibrium gel and an ar- 
rested state generated through spinodal decomposition. 
In the present model, where phase separation is confined 
to the low region of the phase diagram, it is possible to 
reach in equilibrium extremely low T, making it possible 
to study reversible gel formation. 

Besides the low-g growth, on cooling S{q) shows a pro- 
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FIG. 5: (a) Evolution of the static structure factor S{q) with 
T for the (f) = 0.20 isochore for n^ax = 3. Below T = 0.125 
the system has reached an almost fully connected state and 
S{q) does not change any longer with T; (b) Evolution of the 
static structure factor S{q) with (f} for the T — 0.125 isotherm 

for riniax = 4. 

gressive structuring of peaks at qa ~ 2tt and multiples 
thereof, signaling the fact that particles progressively be- 
come more and more correlated through bond formation. 
Indeed, the potential energy of the system progressively 
approaches the ground state value, where all particles 
have nmax bonds js^- Fig. \^ shows the evolution of 
S{q) on increasing (j) along a low T isotherm. Moving 
further from the spinodal, the q ^ peak decreases. A 
small pre-peak, around qcr « 3 is present at low densities, 
and persists with smaller intensity also at intermediate 
(j). Beyond (p ~ 0.45, a significant growth of the nearest- 
neighbour peak is found, signaling the increasing role of 
packing. 

D. Mean squared displacement and caging 

One of the hallmarks of glassy dynamics is the caging 
of a particle by its immediate neighbors. Caging is most 
easily seen in a log-log plot of the the mean squared 



displacement (MSD) vs time as an intermediate-time 
plateau, separating short-time ballistic intracage motion 
and long-time diffusion out of the cage. The height of the 
plateau in the MSD provides a typical (squared) value for 
the localization length Iq of the particles within the ar- 
rested state. For the standard HS glass, Iq is found to be 
roughly O.lcr, and corresponds to the average distance a 
particle can explore rattling within its nearest neighbour 
cage. For an attractive glass, on the other hand, Iq cor- 
responds to the attractive well-width, since particles are 
forced to rattle within the bond distance. For this reason, 
the MSD plateau is significantly smaller than for the HS 
glass (of the order « 10~^ versus (O.Ict)^ « 10~^). 
For the same reason, the q-width of fq is significantly 
larger for the attractive glass solution than for the repul- 
sive one. 

In the low-(/> study reported in 33], on isochoric cool- 
ing a clear plateau develops in the MSD, but its value 
indicates a very high localization length, of the order of 
one or more particle diameters. The localization length 
does not change appreciably with T, even below Tc- We 
attributed this finding to the presence of a long-living 
percolating network of particles allowing for ample sin- 
gle particle movements (arising from the gel "vibrational" 
modes) which completely mask the bond localization. In 
this respect, the connectivity of the network plays an im- 
portant role in the slowing down and provides an addi- 
tional mechanism of arrest with respect to both attractive 
and repulsive glasses. 

Here we study the high-(/) behavior with the aim of lo- 
cating the cross-over from low-^ arrest (gel) to the high-0 
case and to see if a crossover or transition emerges to one 
of the above cited glasses. Results for the (p dependence 
of the MSD along a low T isotherm are presented in Fig.|5| 
for riiiiax = 3 and 4. Both graphs show similar features. 
A very high plateau of order unity, slightly decreasing 
with (j), is found up to (/> w 0.45. Although the long-time 
dynamics is monotonically slower with increasing (j), the 
plateau becomes less defined near (p — 0.50, slowly cross- 
ing over to a quite distinct plateau compatible with the 
HS one at (p — 0.55 for n^ax = 4 and at </> = 0.56 for 

^max — 3. 

Defining the cage length Iq as the square root of the 
MSD value at the inflection point of the MSD in log-log 
scale we plot in the insets of Fig. The cage length 
starts from values larger (n,„ax = 3) or close (nmax = 4) 
to a and progressively approach the HS limiting value 
0.1a. 

In a window of (p values, the crossover between the 
two plateaux in the MSD displays a sub-diffusive behav- 
ior for up to two decades in time, i.e. (r^) cx i", with 
a state-point dependent exponent a < 1. A similar be- 
havior was found in the simulations of the SW system 
[2^ Isij , in a limited T- window, within the liquid reen- 
trant region. In the SW case the subdiffusive behavior 
is found for MSD values between (the bond cage) 
and 10~^a^ (the HS cage) and it is due to a competition 
between attractive and HS glasses at low and high T. Ex- 
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FIG. 6: Mean squared displacement (MSD) and caging, (a) 
MSD along T = 0.1 for n^ax = 3; (b) MSD along T = 0.125 
for nmax ~ 4. The insets show the localization length lo as a 
function of (j>, with dotted line as guide to the eye. Dashed 
lines highlight subdiffusive behavior at intermediate (p. 



plicit MCT predictions have confirmed this feature j^, 
connecting it with the presence of a near-by higher order 
singularity. In the present case it appears that the sub- 
diffusive behavior arises from the competition between 
the very different localization lengths of the gel-arrested 
network and the HS glass. The MSD phenomenology is 
reminiscent of that found in presence of a higher order 
MCT singularity 'so'l. It is possible that the subdiffusive 
behavior (and other features such as the logarithmic de- 
cay of the density auto-correlation functions discussed 
later on) arise generically from the competition between 
two disordered arrested states. 



E. Density Relaxation, Bond Relaxation and 
Non-ergodicity 



We now focus on the behavior of the density autocor- 
relation functions Fq{t). Ref. [s^ called attention on the 
different time dependence of the low and high q window. 
At small q, dynamics slow down significantly and become 



non-ergodic, while at larger q (already on the scale of 
nearest neighbours) dynamics remain ergodic to within 
numerical accuracy. At low density, the non-ergodicity 
parameter for the gel is different from that of either the 
repulsive or attractive glass. We now investigate the ef- 
fect of density on these findings. 

To provide a picture of the behavior of the dynamics 
at low T as a function of (j) we plot in Fig. [7| the density 
correlators for three different values of q, corresponding 
to distances respectively much larger, larger, and com- 
parable to the nearest neighbor distance a. In addition, 
we report the behaviour of the bond correlation function 
<j>B{t) at the same low T, for small and large (p. 

The 0-evolution of the shape of the correlation function 
is particularly complex. At very small g, {qa 1), all 
correlation functions show a clear plateau, followed by 
the a-relaxation process. Interestingly, the plateau value 
has a non-monotonic behavior with 0. It starts from 
an high value at the lowest 4> and decreases down to less 
than 0.1 before increasing again on approaching the hard- 
sphere glass. At the present time, we have no explanation 
for these trends. 

Even more complicated is the gtr ~ 3 case. The plateau 
value first increases with 0, then a reversal of the trend 
is observed at (p — 0.40 where the plateau height starts 
to decrease. Such a decrease persists up to (j) — 0.50, 
after which a distinguishable plateau almost disappears. 
Correlators are higher at comparable times and become 
almost logarithmic for up to three time decades. At = 
0.56, a clear repulsive glassy behavior is recovered. 

Much simpler is the interpretation of the last case, 
qa ^ 7 (and larger q). Here, the standard scenario for 
the repulsive glass is observed, despite the presence of a 
connected long-living network of bonds. The absence of 
any detectable (within our numerical precision) plateau 
at small length scales and low (j) confirms the ability of 
the particles to explore distances smaller than a without 
any constraint. This is an effect of the loose character of 
the network, of the small overall cj) (as compared to the 
typical HS glass values) and of the small local degree of 
connectivity. 

Focusing on the behaviour of 0b (t), we find that the 
curves follow closely, at all studied (j), a simple exponen- 
tial law, i.e. a stretched exponential fit gives an expo- 
nent /3 always close to 1 (tending to 1 with decreasing 
T). From the figures, it is evident that bond relaxation 
is always much slower as compared to density relaxation 
even for very small q. This suggests that, up to at least 
a time of order 10, the density relaxation is coupled to 
the movements of a permanent network which, without 
breaking most of its bonds, is capable of spanning a large 
part of the simulation box. At longer times, the break- 
ing of the bonds enters into play, producing a secondary 
very slow relaxation in Fq (t) , accompanied by very small 
plateau, as it is evident in panels (a) and (b) for small q. 

On increasing n,„ax, the same features for Fq{t) are ob- 
served at a progressively larger q. An example is reported 
in Figl^lfor the case qa ^ 8. Here, around = 0.52 a log- 




FIG. 7: Density autocorrelation functions for T = 0.1 and 
rimax = 3 at various 4> and qa ~ 1,3,7. Also, bond correla- 
tion functions (j>B{t) (dashes lines) at small and large (p are 
reported in all three panels for comparison. 



arithmic decay for about two time decades is observed. 
At this wave-vector, corresponding to lengths smaller 
than the nearest neighbour one, the non-monotonic be- 
havior of the plateau is not present anymore. 

To better grasp the logarithmic behavior observed in 
Figures. [7| and |S1 in Fig. El we show the g-dependence 
of Fq{t) at (j) — 0.54, i.e. the </> showing the most en- 
hanced log{t) dependence, for nmax = 3. We note that 
the best (long- lasting) log{t) dependence is seen in a fi- 
nite window of q-values, roughly between 3 < < 6, 
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FIG. 8: Density autocorrelation functions for T — 0.125 and 
nir 




FIG. 9: Density autocorrelation functions for T = 0.1 , (j> ~ 
0.54 and rimax = 3 at various qa. A log{t) behavior is observed 
for 3 < gcr < 6 for up to three time decades. 



and it covers about three orders of magnitude in time. 
Again, on increasing n,„ax this q-window shifts to larger 
q (approximately between 8 < < 10). 

A possible quantification of the characteristic time 
of the dynamics and of the non-ergodicity parame- 
ter is provided by stretched exponential fits {Fq{t) = 
/q cxp — (i/rq)'^".) of the long time decay. This fit al- 
lows us to extract information on the behavior of the 
non-ergodicity parameter fq, as well as on the stretching 
exponent (3q and an estimation of the relaxation time Tq. 
Still, in a certain region of (p and q values, the decay of 
the correlators is clearly different from a stretched ex- 
ponential, being Fq(t) essentialy linear in log{t). Under 
these conditions we cannot fit the density correlators and 
estimate the non-ergodicity parameter. 

Fig.llUlshows fq along the low T isotherm for rimax = 3 
and 4 for several (p. In both different non-crgodic 

behavior at low and high (p is evident. At low cj), fq is 
largest at g ^ 0, then decays rapidly to zero within a 
range of about 5 in units of qa. With increasing 0, the 
overall height of fq decreases, but a small peak starts to 
form at larger g, which is still of the order of a few qa. 
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FIG. 10: Non-ergodicity factor, obtained from streched expo- 
nential fits, at various 4> for T — 0.1 and rimax = 3 (a) and 
T = 0.125 and n^ax = 4 (b). 



This behavior roughly follows that of S{q), for which the 
low q increase turns to a small peak at finite q (see for 
example in Fig|S|3). At large </> the shape of fq closely 
resembles that of a hard sphere system, with the first 
peak around qa sa 2tt, i.e. the nearest neighbour length 
scale. 

At intermediate 4>, we observe a slightly different be- 
havior between the two studied values of nmax- As dis- 
cussed before, for rimax = 3 and 0.45 < (j) < 0.55, we 
cannot estimate fq from the fits. The reason for this 
is that Fq{t) displays unusual features, like a logarith- 
mic decay at certain q followed by a secondary relaxation 
reminiscent of the gel type. Fig. ^| shows that a sharp 
transition in the q dependence of fq takes place between 
0.45 < 4> < 0.55, which we associate to the gel-to-glass 
transition. For rimax = 4 there seems to be a smoother 
transition although a similar non-monotonic behavior of 
fq with (j) at low q is observed. 

We further note that the relaxation time Tq, extracted 
from the fits or otherwise, monotonically increases with 
(j) despite the non-monotonicity in the plateau. The be- 



FIG. 11: (^-dependence of the density relaxation time Tq at 
qa ~ 1,3,7 versus bond relaxation time tb (triangles), at 
fixed T — 0.1 for nmax=3- Tq are either extracted from 
stretched exponential fits or via the relation Fq(Tq) — 0.01 
{qa ~ 3) and Fqijq) = 0.1 {qa ~ 7). In the latter case, the 
dashed line is a fit of the data to a power law {(j) — (j>c)~' . 



havior of Tq with 



at low T is shown in Fig^] at the 
q- values discussed above for n,nax — 3, together with the 
corresponding behavior for the bond relaxation time tb ■ 
Apart from the fits, Tq can be also conventionally de- 
fined as the time at which the normalized correlators are 
equal to an arbitrary (low) value, which is chosen for 
convenience as 0.1. Doing so, for qa ^ 7, we find that 
Tq follows a typical, for glass-forming systems, power-law 
dependence on 0, also shown in the figure, with critical 
(f>c 0.562 and exponent 7 ~ 2.5. This result for (j>c is 
quite consistent with that extracted from the diffusivity 
fits discussed earlier. However, on lowering g, the situ- 
ation becomes more complicated. At qa ~ 3, the con- 
ventional arbitrary value 0.1 is too high, since the corre- 
lators displays smaller plateau values. Thus, in Fig. 1111 
we report Tq obtained both from the fits and by choos- 
ing a conventional value 0.01. The results are parallel 
to each other, and do not show a power law behavior. 
Rather there seems to be a crossover regime at interme- 
diate (f). However, the relaxation time is monotonically 
increasing at this length scale of observation, despite the 
non-monotonic behavior of the plateau. At even lower q, 
e.g. gcr ^ 1, we cannot define a satisfying finite value for 
Fq{t), below the plateaux at all studied and numeri- 
cally detectable. Also the fits cannot be relied on at high 
densities. As shown in the inset, it seems that Tq does 
not vary strongly with 0, and indeed all correlators (see 
Fig. E^) seem to meet at Fq (t) = at around the same 
value of t, up to (p = 0.52. 

A sensible estimate of tb is simply obtained by 
stretched exponential fits of 4)B(t). tb is always longer 
than the density relaxation time at all q, as remarked 
above. Moreover, the increase of tb upon (/), is rather 
small, indicating that bonds are slightly stabilized by 
crowding, tb is completely decoupled from Tq at large 
q, while it becomes coupled to Tq at small q, the data be- 
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ing almost parallel, and, it seems that in the limit g — > 0, 
the two relaxation times would coincide. Once again, we 
associate the large tb to the very slow secondary relax- 
ation observed for density correlators at small q (see for 
example gcr ~ 3 and 0.50 < (p < 0.55 in Fig. Uja). 

A final observation concerns the behavior of Pq ex- 
tracted from the fits, outside the logarithmic regimes. 
We find at high q {qa > 7) and low 0, that values of (3q 
are larger than 1, actually close to 1.5, associated more 
to a compressed than a stretched exponential pMl90ll9l| . 
However, this value decreases below 1 for high cf) and large 
q, while it tends to be 1 for low (j) and low q {qa < 7). We 
note that care has to be taken when exponents greater 
than 1 are found in Newtonian dynamics, since they could 
arise from undamped motion of clusters (possibly of ro- 
tational origin), or elastic motion within the percolating 
structure. A comparison with a Brownian dynamics sim- 
ulation may help to clarify this issue. 

Additionally, we focus on the T dependence of the den- 
sity correlators. In Ref. 33], we discussed the T depen- 
dence along the lowest isochore (j> = 0.20. There, we 
found that, only at low enough q, there was the emer- 
gence of a gel plateau at low T and that, in this region of 
wavevectors, the density relaxation time Tq as well as the 
bond relaxation time tb both followed an Arrhenius law. 
On increasing density, there is the gel to glass crossover, 
which is clearly visible by looking at the T behavior of 
Fq(t), as shown in Fig. ^1 Up to (j) — 0.45 a clear gel 
plateau is approached. Beyond this value, the anoma- 
lously slow logarithmic decay is observed, at its best for 
(j) = 0.54, reported in figure Eta). Above this (p, e.g. 
(j) = 0.55 in ll2r b). a crossover from logarithmic to stan- 
dard glass regime (i.e., typical two-step relaxation) is ob- 
served. A clear difference between these two graphs is 
the fact that at ^ = 0.54 pure logarithmic decay is ob- 
served up to the lowest T, while at (j) — 0.55 a kink, 
evidence of the nearby glass transition, is present. We 
also note that, interestingly, the shape of the correlators 
in Fig. I12r b) is very reminiscent of Fig. 11 in R.ef . [l^ 
and Fig. 6 in Ref. "251, respectively MCT predictions and 
simulation results for the SW model, for fixed tempera- 
ture in the reentrant region and varying density. Here, 
the analogous plot is reported, at a density within the 
gel to glass crossover, and varying T. Again, this sug- 
gests a close similarity of the features at the gel to glass 
crossover with respect to the SW crossover from repulsive 
to attractive glass. 

Finally, both the relaxation time Tq and the bond re- 
laxation time Tb are found to obey Arrhenius dynamics 
in T at all (j). 



IV. DISCUSSION 

The interpretation of these results is not completely 
straightforward. However, they strongly suggest that the 
system gels at low and forms a glass at high ones. It is 
important to understand what happens in between these 
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FIG. 12: Density autocorrelation functions for rimax = 3 at 
various studied T for (a) (j) = 0.54 and qa ~ 4.4 (maximally 
enhanhced log(t) behavior) and (b) = 0.55 and ga ~ 3 
(interference between log{t) behavior and standard glass-like 
a-relaxation) . 



two regimes, where a new type of relaxation takes place, 
as a result of the competition of the two effects. Simi- 
larly to the simple SW case, where anomalous dynamics 
arises in the reentrant liquid region from the competi- 
tion between attractive and repulsive glass, here again 
anomalous dynamics is generated from the presence of 
two distinct arrested states. However, in the attractive- 
vs-repulsive glass scenario, temperature is the control pa- 
rameter generating a liquid pocket in between. Here 4> 
is the control parameter. We have to bear in mind that 
at T = 0.1 the energy per particlce of the system for 
?^inax — 3 varies from about —1.48 at = 0.20 to about 
— 1.495 at = 0.56, indicating that most of the parti- 
cles (~ 99%) are fully bonded (i.e. they have reached 
the Tiniax limit already). With increasing cj), more neigh- 
bors surround each particle but only nmax of them in- 
teract via an attractive well, the others probing only the 
hard core interaction. Thus we cannot expect a non- 
monotonic (j) dependence of the characteristic time (i.e. 
no re-entrance). 
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Where does the logarithm/subdifFusivity come from? 
A more intuitive understanding of the anomalous dynam- 
ics results from interpreting the MSD behavior. If one 
thinks simply of a filling up of space, the MSD plateau 
should monotonically decrease and no subdiffusive be- 
havior should be observed. However, at « 0.54 a clear 

law is observed. After the ballistic regime, particles 
start to feel the presence of the nearest-neighbours and 
the MSD slows down. At long times, particles are able to 
break and reform the bonds, the network fully restruc- 
tures itself and proper diffusion is observed. In the inter- 
mediate time window one observes the competition be- 
tween excluded-volume confinement and exploration of 
space associated to the motion of the unbroken network. 
We believe that this is at the origin of the anomalously 
slow diffusion and logarithmic decay. 

To support this hypothesis we note that the logarith- 
mic decay shows up only in a window of small q's sig- 
nificantly less than the nearest-neighbour inverse length 
(i.e. over distances where connectivity is probed). On 
increasing rimax, this length becomes smaller due to the 
higher degree of constraint for the network, in the same 
way as the localization length, estimaed from the MSD, 
decreases. However, we recall that in the SW system 
such anomalous log{t) decay was observed for very large 
q, associated to the typical distance of the short-range 
attraction. This provides further evidence that in the 
f^max case the connectivity of the network is associated 
with generating unusual log{t) features and confirming 
that the width of the attractive well does not play any 
significant role. 

A final comment regards the existence of the so-called 
higher order MCT singularity. In the SW case, the width 
of the attraction A is the crucial control parameter driv- 
ing the system close to the singularity. In the nmax case, 
A does not play a relevant role. It is intriguing to ask 
ourselves which parameter plays the role of A, if the gel- 
to-glass cross over belongs to the same class of models 
possessing a higher order MCT singularity. One possible 
answer is nmax itself. Indeed, we noted how the char- 
acteristic length scale to observe logarithmic behavior in 
Fq{t) shifts with increasing n,„ax- Our present knowl- 
edge is that this is found at qcr « 4 for rimax = 3 and 
qa fa 8 for nmax = 4. We also know that for the SW, 
although being associated to the competition of repul- 
sive and attractive glasses, it moves up to qa ~ 20|81|. 
The possible existence of a smooth crossover between the 
two phenomena could be theoretically investigated with 
simulations. However, the existence of the higher order 
singularity in the present model is destined to be uncer- 
tain, unless some theory is devised for the gel transition 
and its predictions tested against the simulations. 



V. CONCLUSIONS 

Understanding the slowing down of the dynamics in 
colloidal systems and the loci of dynamic arrest in the 



full {(j)—T) plane, encompassing gel and glass transitions, 
is one of the open issues in soft condensed matter. Two 
classes of potentials have been explored in some details 
in the recent years: (i) short-range attractive spherical 
potentials and (ii) short-range attractive spherical po- 
tentials complemented with a repulsive shoulder. In the 
first case, it has been shown that low arrested states 
arise only as a result of a interrupted phase separation. 
The second case appears to be much more complicated 
and partially unresolved. For certain values of the pa- 
rameters of the repulsive potential the system separates 
into clusters (which can be interpreted as a mesoscopic 
interrupted phase separation) and dynamic arrest at low 
(f) can follow from a cluster glass transition or from cluster 
percolation. 

In the attempt to provide an accurate picture of dy- 
namics and a model for dynamic arrest at low (j) in the 
absence of phase separation (both at the macroscopic and 
mesoscopic level) we present here a study of a minimal 
model of gel-forming systems. The model builds up on 
the intuition that phase separation is suppressed when 
the number of interacting neighbors becomes less than 
six, since the energetic driving force for phase separation 
becomes less effective jo^- To retain the spherical aspect 
of the potential, a standard SW interaction potential is 
complemented by a constraint on the maximum number 
JT-max of bonded neighbors (a model similar to the one 
first introduced by Speedy and Debenedetti [gI. l65l| ) . 

In this manuscript we have presented a detailed study 
in the (0, T) plane of the dynamics for nmax = 3 and 
f^max — 4. For these two values of n^ax the region of 
phase diagram where unstable states (with respect to 
phase separation) are present shrinks to T < 0.1 and 
(j) < 0.3, making it possible to approach on cooling low 
arrested states, technically in metastable (with respect to 
crystallization) equilibrium. The simplicity of the model 
makes it possible to study it numerically even at very 
low T and estimate, with accuracy, the low-T fate of the 
supercooled liquid. It becomes possible to predict the 
regions in the phase diagram where disordered arrested 
states are kinetically stabilized as compared to the lowest 
free energy crystalline states. 

Dynamics in the rimax model is also important be- 
cause it provides a zero-th order reference system for the 
dynamics of particles interacting via directional poten- 
tials. These systems include globular protein solutions 
(hydrophilic/hydrophobic patches on the surfaces of pro- 
teins), the new generation of patchy colloids, and, at a 
smaller scale, network forming liquids. In this respect, 
the nmax model allows us to study the generic features 
(since it neglects the geometric correlations induced by 
directional forces) of particle association. It has the po- 
tential to provide us with an important reference frame to 
understand dynamical arrest in network-forming liquids 
and the dependence of the general dynamic and thermo- 
dynamic features on the number of patchy interactions. 

One of the important results is contained in Fig. ^ 
which shows the locus of iso-diffusivity in the (0 — T) 
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plane. These lines, which provide an accurate estimate 
of the glass transition line {D 0) are found to be es- 
sentially vertical at high (p, in correspondence to the HS 
glass transition, and essentially horizontal at low T. Only 
a very weak, almost negiiglible, reentrance in (jj is ob- 
served. Extrapolating the D-dependence by power laws 
in (p and by Arrhenius laws in T, we estimate the glass- 
lines. The Arrhenius law is found to be valid for many 
decades, suggesting that, at intermediate (j), D vanishes 
only at T = . Data strongly support the possibility 
of two distinct arrest transitions: a glass of the HS typ e 
driven only by packing at high (f> and a gel at low T [93j . 

We have chosen the word gel to label arrest at small 
and intermediate </> since the analysis of the simulation 
data confirms that the establishment of a network of 
bonded particles and the network connectivity plays a 
significant role in the arrest process. While particles are 
locally caged by SW bonds with nmax neighbors (and in 
this respect one would be tempted to name it an attrac- 
tive glass), particle localization is not only controlled by 
bonding. Bonded particles are free to explore space (re- 
taining their connectivity) until they are limited by the 
network constraints. Indeed, the plateau of the MSD is, 
especially at the lowest 4>, larger than the particle size. 
We remark that the bond localization, typical of the at- 
tractive glass case, is not observed throughout the phase 
diagram, neither in the MSD nor in the width of fq. We 
believe this is due to the fact that, although bonding is 
present, particles are confined by the potential well only 
relative to each other. The network connectivity length 
is the quantity that enters into the determination of the 
localization length in the arrested state. On increasing 0, 
the localization length progressively approaches the one 
characteristic of the hard-sphere glass, signaling that a 
cross-over to the excluded volume case takes place. 

The intersection between the repulsive glass and gel 
loci appears to be associated to anomalous dynamics. 
No intermediate liquid state, i.e. no reentrant regime in 
(f), is found. Interestingly enough, these anomalies are 
strongly reminiscent of the anomalies observed close to 
the intersection of the attractive and repulsive glasses in 
the case of short-range interacting particles. Correlation 
functions show a clear log(t) dependence in a window of 
q vectors and the MSD shows a clear subdiffusive behav- 
ior ~ t". Here, the gel localization length is larger that 
of the HS glass, a different scenario from the attractive- 



repulsive case. These results support the hypothesis that 
a possible MCT-type higher order singularity in the rimax 
model is present and, at the same time, provide fur- 
ther support to the intrinsic difference in the localization 
mechanisms that are active for the two arrested states. 
In contrast to the SW case, the well-width is not a crucial 
length scale in the problem, while an important param- 
eter is n,„ax, that could be the control parameter of a 
putative higher order singularity of the MCT type. 

One further consideration refers to the role of A in the 
"max model. We have chosen to use e = A/(cr -I- A) = 
0.03 to connect with the well studied corresponding SW 
case. We do not expect significant differences for the 
two arrested states for A values up to e ~ 0.1 — 0.2, since 
the connectivity properties would be essentially identical. 
Thus, we still expect the existence of distinct gel and 
glass lines. Only the interplay between the two could 
be affected, as for example the log{t) behavior should 
be shifted in its ^-dependence. Also, in that case, the 
behavior with rimax is not a priori clear, since for larger 
A values no MCT singularity is present for the SW model 
and, from a theoretical point of view, the attractive glass 
line is not physically distinct from the repulsive one. 

In summary, the present model provides a clear indi- 
cation that even if liquid-gas phase separation can be 
avoided and arrest at low (j) can be explored in equilib- 
rium conditions, the observed arrested state is not the 
low-(/) extension of the attractive glass. The present re- 
sults strongly suggest that the attractive glass is an ar- 
rested state of matter which can be observed in short- 
range attractive potentials only at relatively high (p, be- 
ing limited by the spinodal curve. When the inter- 
particle potential favors a limited valency, arrest at low (p 
becomes possible but with a mechanism based on the con- 
nectivity properties of a stable particle network, clearly 
different from what would be the extension of the (at- 
tractive) glass line. 
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